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O , Abstract 

Q ! Analytical formula for multiplicity distribution is derived in the QO ap- 

^ \ proach, where chaotic and coherent fields are contained. Observed charged 

multiplicity distributions in Au+Au collisions at y/s = 200 AGeV and in 
pp collisions at y/s = 900 GeV are analyzed by the formula. Chaoticity 
^ ■ parameters in the inclusive events estimated from the analysis of multiplic- 

ity distributions are compared with those estimated from the analysis of 



p. 



observed two-particle inclusive identical particle correlations. 



^ '. 1 Introduction 

(N . In high energy nucleus-nucleus (AA) collisions or hadron-hadron collisions, 

^ . Bose-Einstein correlations of identical particles are considered as one of the 

^ ! possible measures for the space-time domain where identical particles are 

produced. As the colliding energy of AA collisions increases, thousands of 
identical particles can be produced in an event. Then, the production do- 
main of those particles can be analyzed precisely event by event, or among 
the events with fixed multiplicity. Up to the present, most of theoretical 
approaches to identical particle correlations at fixed multiplicity are investi- 
5^ I gated in the case of purely chaotic field PQ El El II]- 

One of the theoretical approaches to the Bose-Einstein correlations is 
made on the analogy of the quantum optics [5], where two types of sources, 
chaotic and coherent are introduced. A diagrammatical method, based on the 
Glauber-Lachs formula f5j, has been proposed [B] to find higher order Bose- 
Einstein correlation (BEG) functions. In Ref.[7], the generating functional 
(GF) for the momentum densities in the inclusive events is derived, and a 
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diagrammatic representation for cumulants is proposed. Identical particle 
correlations at fixed multiplicity are formulated in [S]. 

In the present paper, analytical formula for multiplicity distribution, 
which is identified to the Glauber-Lachs formula, is derived. A relation be- 
tween the chaoticity parameter psm in the semi-inclusive events and that pi^ 
in the inclusive events is also obtained. Observed multiplicity distributions 
are analyzed by the Glauber-Lachs formula. Estimated value of from the 
observed multiplicity distributions is compared with that from the observed 
inclusive identical two-particle correlations. 

2 Generating functional 

In the QO approach, the n-particle momentum density in the semi-inclusive 
events is defined by, 

PniPl,--- ,Pn) = Co(|/(pi)r--|/(Pn)l'>,, (1) 

M 

f{p) = 5^a,0,(p) + /,(p), (2) 

i=l 

where cq is a normalization factor. In Eq.([2]), (j)i{p) and fdp) are amplitudes 
of the i-th chaotic source and a coherent source, respectively, and is a 
random complex number attached to the i-th chaotic source. The number 
M of independent chaotic sources is regarded to be infinite. 

In Eq.([T]), parenthesis {F)a denotes an average of F over the random 
number with a Gaussian weight [5]: 

M I |2 

(^>" = n(^]^/-p|-^K'>-)'-- (3) 

The single-particle and the two-particle momentum densities are respectively 
given as 

Pi(Pi) = Co(|/(pi)Ha = Co[r(pi,pi) + c(pi,pi)], 

P2{PUP2) = Co(|/(pi)/(p2)Ha 

= co{p{pi)p{p2) + |r(pi,p2)r + 2Re[r(pi,p2)c(p2,Pi)]}, 

where r{pi,p2) is a correlation caused by the chaotic sources and c{pi,p2) is 
a correlation by the coherent source, 

M 

r{pi,P2) = X^Ai0i(pi)0*(p2), C{pi,p2) = fc{pi)fc{P2)- (4) 
i=l 
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The GF of the semi-inclusive events is defined by the following equation, 



Zsm[hip)] = ^ ^ f ■■■ I Pn(Pl, • • ■,Pn)h{pi) ■ ■ ■ /l(p„)^^ ■ ■ ■ (5) 
n=l ^' ^ ^ 



El En 



From Eqs.([T]) and (|5]), the GF is written as 



Zsm[h{p)] = co(exp / I f{p) |2 h{p) 



^ d?p 



(6) 



where an additional constant Zsra[h{p) = 0] is added to the right hand side 
of Eq.([6]). Inversely, the n-particle momentum density in the semi-inclusive 
events is given from the GF as 



PniPl, ■ ■ ■ ,Pn) = Ei---E„ 



S-'Z,Jh{p)] 



6h{pi) ■ ■■6h{pn) 



h(p)=0 



The generating function for multiplicity distribution P{n) is given from 
Eq.([5]) if function h{p) is independent from momentum p: 



,2 d^P-. 



Zsm{h) =co(^expU I f{p) I —h 



The multiplicity distribution (MD) is given from Eq.(I7]) as 

1 d^Z,^{h) 



(7) 



P(0) = Zs^(O) = Co, P{n) 



h=0 



3 Cumulants 



The GF Gsm[h{p)] for cumulants in the semi-inclusive events is defined by 
the following equation. 



Gsm[h{p)] = \nZsm[h{p)]. 
The n-th order cumulant is given by 

S'^g^AKp)] 



gn{pi, . . . ,Pn) = El - ■ ■ Er, 



Sh{pi) ■ ■■6h{pr, 



(9) 



h(p)=0 
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From Eqs.(in]), (E]) and (JH]), we have an iteration relation for momentum den- 
sities, 

Pi(pi) = cogiipi), 

Pnipu...,Pn) = 5'l(Pl)Pn-l(P2, • • • 

n-2 

, Pj-^ , . . . ,pj- 

i=l 

+ Cogn{pu. . . ,Pn)- (10) 

The second summation on the right hand side of Eq. ffTOi) indicates that 
all possible combinations of (ji, ■ ■ ■ , j,) and (ji+i, . . . , jn-i) are taken from 
(2, 3, ... , n). The ra-particle momentum density Pn{pi, ■ ■ ■ ,Pn) {n = 1, 2, . . .) 
is written by the cumulant gi{pi^ ■ ■ ■ ,Pi) (i = 1, 2, . . . , n) from Eq. flTU]) . 

Diagrammatic representation of Eq. (llOI) is shown in [S]. For example, 
cumulants up to the third order are written explicitly as 

gi2 = ri2r2i + Ci2r2i + ri2C2i, 

fi'123 = ^12^^23^31 + Ci2r23r31 + ri2C23r31 + ri2r23C31 + C.C., (11) 

where the following abbreviations are used, 

9ji- -jm = 9m (Pn , • • • , Pj™ ) , = r{pi, pj ) , dj = qpi, pj ) , 

and complex conjugates to the terms explicitly shown are denoted by c.c. 

4 Multiplicity distribution 

From Eqs. fllOj) and fill I) , we obtain a recurrence equation for MD, 

P(n) = -^[Af +jAf\]P(n-j), n = l,2,..., (12) 
i=i 

where, with Ro{ki, /C2) = uji6^{ki — ^2), 

R,{k,k) — , R,ipuP2)= / r(pi,A;)i?,„i(fc,p2) — , (13) 

/d^k r d^k 
Sj-iik,k) , Sj-i{pi,p2)= / c{pi,k)Rj^i{k,p2) . 
CO J CO 

(14) 



4 



In the foUowings, variables are changed from {piLiPiT) to {yi,PiT), where 
Hi is the rapidity of particle i. Correlations r{pi,p2) and c(pi,p2) are both 
assumed to be real, and parametrized as [8], 



riyi,PiT;y2,P2T) = Psm V'pd/i , Pit )p(l/2, P2t) I{^y, ^Pit), 

c{yi,PlT;y2,P2T) = (1 -Psm)A/p(yi,PlT)p(y2,P2T), 



2 /3 „2 

- *-> V I 1 I f 1/ I 

7r3/2 

\2 /A„ \2l 



/(Ay, Apr) = exp[-7L(A?/) - 7t(Apt) 

where Ay = y2 — Vi, ^Pt = P2T — Pit, and Psm is the chaoticity parameter 
in the semi- inclusive events. It is assumed to be constant at present. 
Then, Rj{pi,p2) in Eqs. ffT^ and f lT^ is given by the following form, 

Rjiyi,PiT,y2,P2T) = Njexp[-Aj{yf + yl) + 2Cjyiy2] 



X 



exp[-Uj{p^T + P2t) + 2VrjPlTP2T], 



and 



A-^P^K), Af .^.,(^)-,^, (15, 



where 



^i = | + 7L, v4,-+i = Ai-7i/(A, +a/2 + 7i), (17) 
Ci = 7i, Cj+i = 7LC'i/(A- + a/2 + 7l), 

[/i = ^ + 7T, f/,+i = t/i-7T/(t^.+/3/2 + 7T), (18) 
m = 7T, Wj+i = ixW^/iUj + /3/2 + 7t), 

/vsm\ 0/^3/2, .+1 ^/^-^^/2T^(f/,+/3/2 + 7T) ' 
Recurrence equations, (|T7j) and (fTSi) can be solved Let Fj be defined by 

F,=A,+i + A, j = 0,l,2,..., (19) 
then, is given by a finite continued fraction 0: 

F,- = P,/g„ J = 0,1,2,..., (20) 

P,- = qoP,-i + qiP,-2, J = 2,3,..., (21) 

= qoQj-i + qiQj-2, j = 2, 3, ■ • • . (22) 
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where, go = a + 27L, qi = -7L^ Pq = go, Qo = 1, A = go^ + giCand Qi = go- 
Therefore, Pj and Qj are given respectively by 

Pj = (r2^+' - ri^+')/(r2 - n), Q, = (ra^'+i - n^+')/{r2 - n), 

where 



+ 2^L - a/o^ + 4a7L _ a + 271, + i/a^ + 4a7z, 



ri = 7: , ^2 



2 ' 2 

Then, Aj and in Eqs. lIT^ and (fTBIl are written respectively as, 

r2 -ri 1 + (ri/ra)^ (ri/rs^/^ 



Similarly, f/^ and Wj are respectively given by 

' 2^' ^l-(tlA2)^ ' l-(tlA2)^ 

where 



_ /3 + 27T - V/?' + 4/37T , _ /3 + 27T + + 4/37t 

From Eqs. ffT^ and f l2^ . following expressions are obtained: 

Af = e{l-(ri/r2p}"Hl-(^iA2p}"'' (25) 

Af\ = Aoi^''{l-{n/T2yY''\l-{t,/hy]-\ (26) 

where 

e = ^-^""^f^ = (1 - v^)(l - v^)\^(no), (27) 



v4o = A/r^^V^(l-tiA2)(l-Psm)(rio). (28) 

The generating function for multiplicity distribution P{n) is given by 

00 

U{z) = Z,^{z) = Y,P{n)z^. 

j=0 

Then the differential equation for Il{z) is obtained from Eq. fll2l) as, 

|n(.) = J](Af + jAf_\ ).-^n(.). (29) 



For j » 1, Eqs.( l25|) and ( 126|) can be approximated as A^^^ ~ ^-^ and 
Aq^^^^. Therefore, with the boundary condition 11(1) = 1, we obtain 



(z-l) 



U{z) = {1 -Pin(n)(^ - 1)}- exp[(l-pi,)(n) ^_^^^^^^^^_^^ 



(30) 



where the average multiphcity (n) and the chaoticity parameter in the 
inclusive events are given by the following equations, 



(n) = e/(l - + ^o/(l - 0', Pin{n) = - 0- 
The multiplicity distribution P{n) is given from Eq. (l30|) as 

1 d''' 



P{n) = -—U{z] 

(1+Pin(n))"+1 



exp 



1 +Pin(^) J "V 1 +Pin(n) 



(31) 
(32) 

(33) 



where Ln{x) denotes the Laguerre polynomial. Equation fl33|) is called the 
Glauber-Lachs formula O |10]. The KNO scaling function of the Glauber- 
Lachs formula is given by. 



1 

— exp 

Pin 



Pin 



loi a/(1 -Pin)^ 



(34) 



where Io{z) is the modified Bessel function. 



5 Analysis of experimental data 

At first, observed multiplicity distributions in Au+Au collisions at ^/s = 200 
AGeV [12] are analyzed by the scaling form 

P{n) = (p{z)/{n), z = n/{n). (35) 

Results for Au+Au collisions and the estimated parameters from the analysis 
are shown in FigiU and Table [H respectively. 

The second order BEC function in the QO approach for the Gaussian 
source function is given by 

^2-/^BG _ ^ ^ 2pi„(l - p,^)E2B + PinE2B\ E^B = CXp [-r^g,^,^] _ 



7 



At Qinv = 0, N"^' /N^'^ — 1 = Pin(2 — Pin). It corresponds to correlation 
strength A. 

The correlation strength A, estimated from the experimental data with 
full Coulomb corrections for 0.2 GeV/c < pr < 2.0 GeV/c for 0-30% cen- 
trahty, is about 0.39 [13], which is much larger than Pin(2 — pin) in Table [1] 
estimated from the multiplicity distributions. 
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Figure 1: Analysis of charged mul- 
tiplicity distributions observed in 
Au+Au collisions [I2] by Eq.([35D 



Figure 2: Analysis of charged mul- 
tiplicity distributions observed in 
pp collisions [H] by Eq.( l35|) 



Table 1: Estimated parameters in the analysis of charged multiplicity distri- 
butions observed in Au+Au collisions at y/s = 200 AGeV |l4j . 



Centrality {n) ± 6{n) ± 5pi 



00 - 05% 62.3 ± 0.04 
20 - 25% 34.0 ± 0.03 
50 - 55% 9.31 ± 0.03 



(9.81 + 0.07) X 10-3 
(2.19 + 0.01) X 10-2 
(7.52 + 0.09) X 10-2 



XiJn.d.f Pin(2-Pin) 

879.3/(58-2) 0.0195 
475.9/(48-2) 0.0433 
914.7/(24-2) 0.145 



Then observed multiplicity distributions in pp collisions at a/s = 900 
GeV [H] are analyzed by Eq. (l35|) . Results and estimated parameters are 
shown in Figj2] and Table [2], respectively. 

The correlation strength A, estimated from the experimental data for 
0.1 GeV/c < Pt < 0.25 GeV/c for \r]\ < 0.8, is 0.628 + 0.133 [15j. It is 
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not inconsistent with Pin(2 — Pin) = 0.729 in Table [2] estimated from the 
multiphcity distributions for pseudorapidity range |?7| < 1.0. 

Table 2: Estimated parameters in the analysis of charged multiplicity distri- 
butions observed in pp collisions at = 900 GeV [H]. 

77 range {n)±6{n) ± XmiJn.d.f Pin(2-pin) 

|r7|< 0.5 3.82 ± 0.05 0.579 ± 0.03 45.2/(24-2) 0.823 

H < 1.0 7.59 ± 0.06 0.479 ±0.01 48.5/(42-2) 0.729 

|r;| < 1.3 9.92 ± 0.08 0.434 ± 0.01 71.8/(48-2) 0.680 



6 Summary and discussions 

The analytical formula for multiplicity distribution in the QO approach 
is derived. It becomes the Glauber-Lachs formula. A relation between 
the chaoticity parameter in the inclusive events and that psm in the 
semi-inclusive events is obtained. Multiplicity distributions observed in the 
Au+Au collisions at ^/s = 200 AGeV and in pp collisions at ^/s = 900 GeV 
are analyzed by the scaling form of Glauber-Lachs formula. The correlation 
strength calculated with pi^ is compared with that measured from the second 
order BEG data. 

In Au+Au collisions, the former calculated with pin is much smaller than 
the latter measured from the BEG data. It would be caused by the fact that 
the centrality region for the data sample of MD is different from that for 
the second order BEG. It would be very interesting, whether both values are 
consistent or not, if the MD and the second order BEG are taken from the 
same centrality region. 

In pp collisions, the former is not inconsistent with the latter. 
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